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Relaxor ferroelectrics have been a focus of intense attention due to their anomalous dielectric 
characteristics, diffuse phase transitions, and strong piezoelectricity. Understanding the structure 
and dynamics of relaxors heis been one of the long-standing challenges in solid-state physics, with the 
current model of polar nanoregions in a non-polar matrix providing only a qualitative description 
of the relaxor phase transitions. In this paper, we investigate the local structure and dynamics 
in 75%PbMgi/3Nb2/303-25%PbTi03 (PMN-PT) using molecular dynamics simulations and the 
dynamic pair distribution function technique. We show for the first time that relaxor transitions can 
be described by local order parameters. We find that structurally, the relaxor phase is characterized 
by the presence of highly anisotropic correlations between the local cation displacements. These 
correlations resemble the hydrogen bond network in water. Our findings contradict the current 
polar nanoregion model; instead, we suggest a new model of a homogeneous random network of 
anisotropically coupled dipoles. 



Recently, relaxor ferroelectrics have become important 
in technological applications, resulting in a revival of in- 
terest in this longstanding fundamental scientific prob- 
lem. [IH5] Compared with normal ferroelectrics, relaxors 
exhibit a stronger piezoelectric effect, a high permittiv- 
ity over a broad temperature range, and unique dielec- 
tric response with strong frequency dispersion. The in- 
verse dielectric response starts deviating from the Curie- 
Weiss law at the Burns temperature (Tf,), significantly 
above the Curie temperature (Tc). For the past several 
decades, these effects have been ascribed to the appear- 
ance of polar nanoregions (PNRs) which form spherical 
or elliptic clusters in non-polar matrix at Tb due to ran- 
dom fields in the material, with the size and the interac- 
tions of the PNR increasing as T is lowered to the Vogel- 
Fulcher freezing temperature (Tf). However, this model 
provides only a qualitative description of the changes in 
the structure through relaxor transitions. Furthermore, 
recent studies using Raman, NMR, neutron scattering 
pair distribution functions (PDFs), and diffuse scatter- 
ing techniques [THU] have demonstrated that static local 
polarization on at least a nanosecond time scale appears 
only at a temperature T* roughly halfway between Tf, 
and Tf. 

In this work, we use molecular dynamics simula- 
tions and analysis of PDFs for 75%Pb(Mgi/3Nb2/3)03- 
25%PbTi03 (PMN-PT) to show that relaxor transitions 
are characterized by well-defined and observable local or- 
der parameters and are due to the onset of anisotropic 
nanoscale correlations of the in-phase cation motions. 
These correlations do not form clusters and therefore can- 
not be explained by the current PNR model. Rather, 
the couplings between displacements are analogous to 
the network of hydrogen bonds in water. Both in wa- 
ter and in relaxors, the presence of such a network in 
a polar environment leads to unique physical properties. 
We therefore propose that the PNR model should be re- 
placed with a model of a hydrogen-bond-like network of 
dipoles generated by local anisotropically coupled cation 




FIG. 1: (Color online) The frequency averaged DPDF. (a- 
b) Color contour plot of experimental G(r, cj) for PMN 
integrated from 10 meV to 20 meV and from -5 meV 
to 5 meV. [10| (c-d) Color contour plot of computational 
0.75PMN-0.25PT G(r, cj) integrated from 10 meV to 20 meV 
and from to 5 meV, respectively. Arrows mark the r values 
for which a ID plot of DPDF intensity versus temperature 
is presented in (e) integrated G(r, w) for the 10-20 meV fre- 
quency range at r=2.38 A, (f) for the 0-5 meV frequency 
range at r=3.3 A, and (g) for the 0-5 meV frequency range 
at r-=2.58 A. See text for interpretation. 



displacements. 

We study an 8640-atom superceU of PMN-PT with 
bond-valence (BV) molecular dynamics (MD) simula- 
tions [H [3j in this work. Application of the dynamic 
pair distribution function (DPDF) method [TD] enables 
identification of the size and directions of displacement 
correlations and the atomic-scale local order parameters 
for relaxor transitions (Computational details of MD and 
methodology of DPDF are described in supplementary 
information) . 

We first compare our results with the experimen- 
tal 0.75PMN-0.25PT phase transition temperatures and 
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with integrated experimental DPDF data of Dmowski 
et al. [10| for PMN(Figure la-b). We integrate the total 
DPDF from 10 meV to 20 meV and from meV to 5 meV 
for 2 < r < 5 A for our 0.75PMN-0.25PT material (Fig- 
ure Ic-d). For the higher frequency interval (Figure Ic), 
a short Pb-0 peak appears at ~550 K, same as the exper- 
imental Tb=550 K of 0.75PMN-0.25PT. 13] This can be 
clearly seen in a ID plot of the integrated high-frequency 
DPDF intensity at r=2.38 A versus temperature (Fig- 
ure le) , where a more rapid rise of DPDF intensity with 
lower T sets in at 550 K. Though the color change in 
Figure Id is subtle, the ID plot for the integrated lower- 
frequency DPDF intensity at r=3.3 A (Figure If) shows 
two changes of slope at 550 K and at 425 K. The tran- 
sition at 425 K is at a temperature that is intermedi- 
ate between the experimental rf,=550 K and T/=380 K 
of 0.75PMN-0.25PT. [17] The integrated lower frequency 
DPDF intensity for the short Pb-0 distance (Figure Id) 
shows a split of the Pb-0 peak away from the main Pb- 
O and 0-0 distances peak at 2.83 A. Since the splitting 
gives rise to a decrease of the DPDF intensity at the dis- 
tance between the two peaks, we examine the integrated 
lower frequency DPDF intensity for r = 2.58 A (Figure 
Ig). This shows two changes of slope at 525 K and at 
350 K, with the rapid intensity decrease and strong peak 
splitting appearing at temperatures slightly below the ex- 
perimental T/=380 K of 0.75PMN-0.25PT. Similar fea- 
tures are present in the integrated DPDF intensities for 
PMN obtained by Dmowski et al. (Figure la-b). Here, in 
integrated DPDF in the higher frequency interval (Figure 
la), the 2.4 A peak appears at sa600 K, close to rb=630 K 
of PMN. In the lower frequency (Figurelb), at r = 3.3 A 
there is a strong intensity enhancement at ^300 K, be- 
tween PMN Tt=630 K and r/=200 K, jTU] and a faint 
increase in intensity above 300 K. The 2.4 A peak visibly 
splits off from the main 2.8 A peak at ~170 K, slightly 
below T/=200 K. The correspondence between the tran- 
sition temperatures found in our DPDF data and the ex- 
perimental 0.75PMN-0.25PT values as well as the good 
agreement between the main features of the experimental 
PMN and computational 0.75PMN-0.25PT DPDFs show 
that our simulations are a good basis for investigating the 
structure and dynamics in relaxor ferroelectrics. 

Cation-oxygen instantaneous and time-averaged PDFs 
(Figure 2) reveal the changes in the dynamics as the 
temperature is lowered and the system undergoes a se- 
quence of transitions. The peak of the instantaneous Pb- 
O partial PDF g{r,t = 0) presented in Figure 2a cor- 
responds to the Pb-0i2 cage and is asymmetric for all 
temperatures with the positively skewed distributions of 
g{r,t — 0) above T = 550 K. This means that Pb atoms 
shift away from the center of their O12 cages even in the 
high-temperature paraelectric phase, in agreement with 
previous results. [H [THl HH] Unlike g{r,t — 0), the time- 
averaged Pb-0 partial PDF G(r, w = 0) in Figure 2b 
shows almost symmetric Pb-0 peaks centered at 2.81 A 
for T >475 K, indicating a small and temperature in- 
dependent time-averaged Pb atom off-centering. Such 





"100 200 300 40n 500 600 700 
Temperature (K) 



100 200 300 400 500 600 700 
Temperature (K) 



FIG. 2: (Color online) Temperature dependence of instanta- 
neous and time averaged PDFs for cation-O atomic pairs, (a) 
Instantaneous g{r,t=0) and (b) time-averaged G{r,uj=0) Pb- 
O partial PDFs for the first Pb-0 peak, (c) The magnitude 
of the instantaneous Pb off-center displacement and its static 
and dynamic components as a function of temperature, (d) 
Nb and Ti displacements as a function of temperature as ob- 
tained from the Nb-O and Ti-0 g{r,t^O) and G(r,a;=0). See 
text for interpretation. 



small time-averaged local cation displacements were ob- 
served in previous simulations and are due to the strong 
local random fields created by the local variations in the 
_B-cation arrangement. p[ IT5] 

To clearly show the changes of local structure with 
temperature, we use the short Pb-0 peak positions for 
g{r,t — 0) and G{r,uj = 0) to calculate the magnitude 
of the instantaneous and the time-averaged static local 
Pb displacements (-Dp^' and Dp^'"^) shown in Figure 2c. 
These are obtained by subtracting the peak position from 
2.83 A Pb-0 distance of the high symmetry Pb-0i2 cage 
at the 0.75PMN-0.25PT lattice constants. We also show 
the dynamic component of the Pb displacement i?p^" 
defined as the difference between Dp^* Dp^'"^. 

Examination of the data in Figure 2c shows that three 
phase transitions occur at rf,=550 K, r*=450 K and 
Tf=375 K, separating four distinct regions. The values 
of all three transition temperatures are agreement with 
experimental 0.75PMN-0.25PT data. P^ITT] For T > Th. 
Dp^* and Dp{^*"^ are both small and change little with 
temperature. For 450 K< T <550 K, Dp^'^ increases 
rapidly as T is lowered, while Z?p^'"^ is the same as for 
T > 550 K. The increasing dynamic displacement com- 
ponent indicates the onset of the dynamic relaxor 
phase. For 400 K< T <450 K, Dp^^^" increase slowly and 
^Pb" starts to plateau. Finally, at r=375 K, another 
transition takes place with Dp^*"^ increasing rapidly and 
D'^ decreasing and both then saturating at their low 
temperature values, as the system undergoes a transition 
into the frozen phase. 

The transition at T* is characterized by the change 
in £)p|i^*"^; we therefore assign the difference between 
the calculated Dp^*''^ and the Dp,^'"^ of the paraelectric 
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phase as the order parameter for the low-teraperature 
phases at T <475 K. The Dp^^''^ is related to the 
Edwards- Anderson spin glass order parameter q, defined 
as q—T,i < Si >^ PD]: which measures the average mag- 
nitude of the static local polarization in the material and 
has been used to model relaxor behavior. [9l[T8l[2T] Using 
NMR experiments on PMN and an analytical random- 
bond random-field spin model, Blinc et al. have pre- 
viously shown that q rises above zero at a temperature 
between Tf, and Tf. [9] This is similar to our finding of 
£)static ]-jgjj^g above the paraelectric value at T*=450 K. 
Since q and Dp^^*"^ do not change at T^, the freezing of the 
local displacements that has been the focus of previous 
investigations is not the relevant process for the onset 
of relaxor behavior at Tf,. Rather, the transition at Tf, 
is characterized by the increase in the local scalar mag- 
nitude of the instantaneous cation displacements D™^ 
(Figure 2c). We therefore assign the difference between 
the calculated and the -Dp^* of the paraelectric 

phase as the order parameter for the dynamic relaxor 
phase found for «475 K< T <550 K. 

Inspections of Nb-0 and Ti-0 partial PDFs (Figure 
2d) show trends similar to those observed for Pb displace- 
ments (Mg cation displacement magnitudes are close to 
zero for all temperatures). The main difference between 
the Pb and i?-cation displacements is that the dynamic 
component of the off-center displacement T)'^^" peaks at 
a higher temperature for the S-cations than for the Pb 
atoms. The cation-oxygen bonding is more distributed 
and flexible in the PbOi2 cuboctahedron than in the BOq 
octahedron. This makes the freezing temperature of the 
Pb displacements lower than for the Nb and Ti cations. 

The fact that a local quantity such as Pb displacement 
magnitude shows order parameter behavior implies that 
local interactions are at the root of relaxor behavior. En- 
hanced local correlation between the Pb displacements is 
the physical origin of the increase in -D"^* for T < 550 K. 
Pb atom displacements in neighboring unit cells are cou- 
pled, so that correlated displacements in the same di- 
rection (even if time-averaged to zero) allow a greater 
magnitude of individual Pb ion off-centering. We now 
examine the cation-cation DPDFs to reveal the changes 
in displacement-displacement correlations that give rise 
to the structural and dynamical properties of the relaxor 
phase. 

We compare Pb-Pb G{r,u}) along the (100), (110) and 
(111) high symmetry directions. In Fig. 3, we show only 
the (100) and (110) directions since the intensities along 
(111) directions are only slightly weaker than the inten- 
sities along (110) directions. The G{r,uj) reveal the spa- 
tial extent and the frequency spectrum of the correlated 
cation motion. We find that the correlations between 
Pb displacements are enhanced between Tf and Tf, and 
exhibit a strong direction dependence. Starting at Tf,, 
the G{r,uj) data along the (100) direction show a strong 
increase in intensity at low w, with a shift of the low 
frequency peak to below 0.1 meV as T approaches Tf. 
The intensity weakens and the peak position changes as 




FIG. 3: (Color online) Pb-Pb G{r,uj) along the high- 
symmetry directions as functions of temperature (K) and fre- 
quency (meV) at the peak position of the instantaneous PDF 
for each temperature. The contour is projected at bottom 
with a color bar which shows the intensity of the peaks. Data 
for (a) r=(100) (b) r=(500) (c)r=(110) (d)r=(330) are shown. 

r increases from 4 A to 20 A, but the appearance of low- 
uj peaks at Tf, and their shift to lower frequencies with 
lower T is present for all DPDF along the (100) direction. 
By contrast, although the (110) and (111) directions also 
show peaks in the in-phase vibration intensity, they de- 
cay dramatically with increasing distance. For example, 
for all temperatures, we find essentially zero DPDF in- 
tensity for the (330) Pb-Pb peaks, which indicates a lack 
of correlated in-phase oscillations. This shows that be- 
low Tf,, the coupling between the local dipoles created by 
Pb displacements is anisotropic, with strong interactions 
only between the dipoles located along Cartesian direc- 
tions. Anisotropic correlations between cation displace- 
ments were also very recently reported by Akbarzadeh et 
al. in Ba(Zr,Ti)03 relaxor. [3T] The direction depen- 
dences of the B-cation-iS-cation correlations are more 
substantial and anisotropically strong coupling, confined 
to B-cation sublattice, along the Cartesian directions is 
present even in the high temperature paraelectric phase 
ofPMN-PT. (See in supplementary information) . This is 
in disagreement with the current model of PNR inside a 
non-polar matrix, where correlations should extend along 
all directions as the temperature is lowered below Tf,. 

The stronger correlation along the (100) direction in- 
dicates that the Pb displacement coupling is not solely 
due to the dipole-dipole interactions. The anisotropy 
is induced by the through- oxygen interactions between 
nearest-neighbor cations that share one (for the id- 
eations) or more (for Pb) O atoms along a Cartesian 
direction in the ABO3 structure. These interactions are 
stronger for the i3-cations, due to the much higher av- 
erage bond valence of each i?-cation-oxygen bond (4/6) 
compared to the valence of the average Pb-0 bond (1/6). 
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The sharing brings strong B-cation displacement cou- 
phng even in the paraelectric phase. On the other hand, 
for Pb atoms, the couphng along (100) is weaker and is 
therefore absent above Ti,. 

Despite the presence of strong i3-cation-i3-cation dis- 
placement correlations, at T > Tf, PMN-PT still exhibits 
normal paraelectric behavior. We ascribe this to the 
fact that the strongly correlated B-cation displacement 
chains are one dimensional and therefore cannot undergo 
a phase transition into a more ordered phase. In ad- 
dition to the through-oxygen interactions, the chains of 
the i?-cation-i?-cation correlations enhance Pb atom cor- 
relations along the Cartesian directions. At T^, as the 
Pb atom displacements next to the strongly correlated 
B-cation chains become correlated, the correlated region 
becomes three-dimensional, enabling a phase transition 
from a disordered paraelectric phase to a locally corre- 
lated dynamic relaxor phase. 

Our results show that the region of correlation is 
strongly anisotropic, with a network of tube-like corre- 
lations along the Cartesian axes and a radius of around 
7.5 A in the other directions. This is the same as ~ 7 A 
estimated for the size of PNR by Gehring et al. based on 
the "waterfall effect" in 0.8PMN-0.20PT. [3] The change 
of the low Lo part of the spectrum from a sharp peak at 
T > Tf, to a broad band absorption for T < Ti, is in 
agreement with the waterfall effect observed by neutron- 
scattering experiments in relaxor materials. A water- 
fall dispersion relation will result in a density of states 
with a peak at higher ui and a flat region on the low- 
frequency side of the peak, precisely as observed in our 
G{r,uj) for T < Ti,. Therefore, our simulations suggest 
that the experimentally observed waterfall effect is due to 
the correlated in-phase vibrations of cation pairs coupled 
by short-range through-oxygen interactions. 

The correlations between cation displacements show a 
strong dependence on the identity of the atoms. For low 
id, high intensity of in-phase vibrational correlation is fa- 
vored by Nb and is disfavored by Mg, with Ti intermedi- 
ate. This is particularly pronounced for the Pb-B-cation 
G{r,uj) where the Pb-Mg G(r, w) shows weaker changes 
with temperature for the short r=3.5 A (111) peak (See 
in supplementary information), in contrast to all other 
atom pairs at r < 4A (including Mg-Mg). Previous 
BVMD simulations have found that Mg atoms and Nb 
atoms surrounded by Mg atoms (Nb^'^) exhibit fast Ar- 
rhenius dynamics typical of the high-temperature para- 
electric phase even for T < T/. [3 [52 Since Mg-Mg, 
Mg-Nb and Mg-Ti G{r, uj) all exhibit the intensity at low 
id associated with relaxor behavior, we assign the lack 
of Pb-Mg coupling as the cause of the paraelectric-like 
dynamics of the Mg and Nb^'^ atoms. 

The current understanding of the relaxor transitions is 
that at Tf, small and dynamic spherical PNR form within 
a paraelectric matrix; as T is lowered, the PNR grow, 
show smaller polarization fluctuations, and start to freeze 
in at T*. [TT1|211I2S] At Tj, a percolation transition takes 
place as all of the PNR coalesce into a a single cluster re- 



sulting in the freezing of the local polarizations. The 
static relaxor state is described as inhomogeneous and 
consisting of static PNR in a statically non-polar matrix, 
with volume fraction of the PNR of about 0.3. [TB] 

Our results cast doubts on the conventional picture of 
relaxor structure and dynamics. Recent MD simulations 
have shown that structural and dielectric properties of 
PMN-PT can be reproduced without PNR. 01 [18; Re- 
laxor behavior in Pb-based perovskites was also shown to 
be controlled by local structure parameters. [57] In this 
work, we find that the average magnitude of local cation 
displacements is quite large at T < r/=:375 K, with Dpb 
magnitude of «0.40 A, close to that found in the pro- 
totypical normal ferroelectric PbTiOa. In ferroelectric 
perovskites, uncorrelated displacements incur a large en- 
ergy cost due to oxygen atom under- and overbonding, 
increased A-B cation repulsion, and unfavorable dipole- 
dipole interactions. [551 [IS] Therefore, a large magnitude 
of local polarization requires strong correlated displace- 
ments, as found in our simulations of PMN-PT; this is 
inconsistent with the idea that most of the material exists 
in a non-polar matrix state. Additionally, we find that 
correlations between the cation displacements are highly 
anisotropic and are weak for (110) and (111) directions 
even at temperatures close to Tf ; this is inconsistent with 
the picture of the strongly polarized spherical nanoclus- 
ters. 

We therefore suggest a alternate model for relaxor 
structure and dynamics. The transition at Th is charac- 
terized by a shift from dipole-dipole interactions to the 
anisotropic short-range, through-oxygen coupling along 
the Cartesian directions. At high T, the thermal energy 
is high enough to disrupt the through-oxygen coupling, 
so that dipole-dipole interactions play a dominant role 
in determining the dynamics of the system. At Tf,, the 
energy lowering due to correlated displacements and the 
decrease in oxygen atom overbonding is larger than the 
entropy cost of correlated displacements. Therefore, the 
structure changes to a network of strong, through-oxygen 
coupled motions along the Cartesian axes. Such a tran- 
sition is similar to the changes that take place when su- 
perheated water is cooled down to room temperature, 
where a standard polar liquid local structure dominated 
by dipole-dipole interactions is transformed into a direc- 
tional, highly anisotropic random H-bond network with 
unique structural and dynamical properties. Extending 
an analogy suggested by Pirc and Blinc relating the fer- 
roelectric phase to solids and the paraelectric phase to 
liquids, [^ we suggest that relaxors are analogous to 
hydrogen-bonded water. A coupling network mediated 
by O atoms plays the role of the H-bonds in water and 
preferential bonding directions lead to dynamic clusters 
of correlated displacements of various sizes. The fact 
that water exhibits Vogel-Fulcher dielectric response in- 
cluding some extremely slow relaxation processes due to 
the collective H-bond network behavior further supports 
this analogy. [3^ 
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Supplementary information 

SI. COMPUTATIONAL DETAILS 

We use a 12x12x12 0.75P(Mg(i/3)Nb2/3)- 
0.25PbTiO3(PMN-PT) supercell with random-site 
ordered B-cation arrangement ^ and perform MD sim- 
ulations at a range of temperatures from 100 K to 700 K 
using an atomistic bond-valence potential [5] derived 
from first-principles calculations and used successfully 
in a previous study of 0.75PMN-0.25PT. [3] The BV 
potential parameters in this work are shown in Table SI 
and S2, with the Og tilt angle potential parameter set 
to 1.6875meV/(deg)^. We run our simulations for 0.5 ns 
in order to be able to extract low frequency data. 





IP (e) 


Sp (eV) 






Pb 


1.574 


0.0618 


2.044 


5.50 


Mg 


0.185 


1.1555 


1.622 


4.29 


Nb 


0.503 


0.3874 


1.907 


5.00 


Ti 


0.575 


0.5839 


1.806 


5.20 


O 


-0.672 


0.4976 


N/A N/A 



Table SI: The charges, qp, and bond-valence model pa- 
rameters, Sfi, for the atomistic model of 0.75PMN-0.25PT. 
The values of r^^ and Cppi we re empirically determined by 
Brown.[4] The notations accord with Reference2. 



Pb Mg Nb Ti O 

Pb 2.221 0.553 2.454 2.312 1.676 

Mg 0.553 1.344 2.920 1.442 1.415 

Nb 2.454 2.920 0.408 3.053 1.278 

Ti 2.312 1.442 3.053 1.037 1.181 

O 1.676 1.415 1.278 1.181 1.927 

Table S 2: The parametrized values of coefficients in the 
interatomic repulsive potentials, Bpfji in Reference 2, for our 
0.75PMN-0.25PT model. 



where N is the number of atoms, 6 is a neutron scatter- 
ing length, (6) is the average of the scattering lengths, a is 
the normalization factor of the Gaussian smearing width, 
and R is an atomic position. The g{r,t) can be Fourier 
transformed to G{r,uj) which denotes DPDF. This makes 
it possible to resolve the vibrational frequencies for the 
different types of atomic pairs at different distances, re- 
lating G(r, u) results to the experimental studies of dy- 
namics. We use three different PDFs in this work for 
structural and dynamics analysis: g{r,t=0), G{r,uj — 0) 
and G{r,uj). The integral of G{r,uj) over w ^ gives 
the difference between g{r,t—0) and G{r,uj = 0). The 
physical interpretation of these different PDF types can 
be illustrated as follows. Consider a generalized PDF 
9ab t) defined as the probability of an atom of species 
A at time t' located at distance r from the position of 
an atom of species B at time t" , such that t' — t"~t. 
Then gABi^jt = 0) is the instantaneous PDF. A peak in 
gAB{r,t = 0) is broadened due to structural variation in 
the preferred interatomic distance r and due to out-of- 
phase vibrations of the A — B pairs around the minimum 
energy r. The in-phase vibrations preserve the preferred 
distance r and do not broaden the peak. The widely used 
elastic neutron scattering technique measures the instan- 
taneous PDF. Similar to the peaks in gABi^^t 0), a 
peak in Gab is broadened by structural variation in the 
preferred interatomic distance r and by out-of-phase vi- 
brations of the A~ B pairs around the minimum energy 
r. Additionally, it is also broadened by the in-phase vi- 
brations oi A — B pairs. This means that a difference 
between gAB{r,t — 0) and GABifjiv = 0) specifies the 
amplitude and the intensity of the in-phase vibrations of 
the ions. The total difference g ~ G can be decomposed 
into contributions at each uj. Thus, DPDF analysis iden- 
tifies at which frequencies particular types of atom pairs 
at particular distances are vibrating in-phase. This ca- 
pability is particularly relevant to the study of relaxor 
dynamics where the appearance of in-phase vibrations of 
the cations characterizes the transition from the para- 
electric to the relaxor phase. 



S2. DYNAMIC PAIR DISTRIBUTION 
FUNCTION ANALYSIS METHODOLOGY 

In this section, we describe procedure of dynamic pair- 
distribution function(DPDF) calculation and explain its 
physical meanings and usefulness. Using MD simulation 
results, we then calculate the g{r,t) generalized PDF 
which measures the probability of the position of one 
atom at time t' being separated by a distance r from the 
position of the other atom at time t" , such that t' — t"—t. 
The (7(r, t) denotes the correlation function at a given dis- 
tance r and time delay t: 

1^1 1 
e dt', (1) 



S3. PB-PB DPDF ALONG THE (111) 
DIRECTION 

The DPDF, G{r,uj), of Pb-Pb atomic pairs for r=(lll) 
and r=(222) are shown in Fig.^ Similar to the DPDFs 
along the (110) direction, there is a rapid decay in peak 
intensity with distance. 



S4. NB-NB DPDF 

DPDF peaks, G{r,uj), for Nb-Nb atomic pairs are 
shown in Fig.^ The G{r,uj) for the (400) Nb-Nb peak 
shows strong positive intensity, whereas G{r, w) for the 
(220) and the (222) Nb-Nb peaks show almost no pos- 
itive intensity despite the shorter interatomic distance. 
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However, here the intensity for the (nOO) peaks is quite 
strong aheady in the paraelectric phase. This is different 
from the behavior observed for the Pb atoms in Figure 3 
of main text, where in the paraelectric phase the in-phase 
vibrations are only weakly enhanced for the (nOO) peaks 
compared to the (nnO) and the (nnn) peaks. 




Fig. SI: (Color online) First and second neighbor Pb-Pb 
DPDF along (111) directions as functions of temperature (K) 
and frequency (meV) at the peak height of the instantaneous 
PDF for each temperature. The contour is projected at bot- 
tom and a color bar the intensity of the peaks. Data for (a) 
r=(lll) (b) r=(222) Correlations of in-phase low-frequency 
vibration along the (110) and the (111) directions at r > 8 A 
get weaker even than the correlations at r=(500). 




Fig. S 2: (Color online) Nb-Nb DPDF along the high- 
symmetry directions as functions of temperature (K) and fre- 
quency (meV). Distance corresponds to the peak position of 
the instantaneous PDF for each temperature. The contour 
is projected at bottom with a color bar which shows the 
intensity of the peaks. Data for (a) r=(200) (b) r=(400) 
(c)r-=(110) (d)r-=(220) are shown. For the (100) direction, 
enhanced intensity at low uj for T < T;, is present even for 
large distances. For the (110) and the (111) directions low- 
frequency vibrations disappear for r > 8 A. In-phase vibra- 
tional intensity at higher uj is more intense for the Cartesian 
peaks than for non-Cartesian peaks for all T. 



S5. PB-MG DPDF 

The G{r,uj) peaks for r=(lll) Pb-Mg atomic pairs are 
shown in Fig.^ Despite a shorter distance, the en- 
hancement of intensity at low lo for T < T;, is weaker 
than for Pb-Pb G{r,uj) at r = 4 A. Such weak coupling 
is unique compared with other nearest-neighbor cation- 
cation atomic pairs in the PMN-PT. 




Fig. S3: (Color online) Pb-Mg DPDF peaks as function of 
temperature (K) and frequency (meV). The intensity values 
taken at the peak position of the first Pb-Mg peak in the 
instantaneous PDF (r ~ 3.5 A) at each temperature. The 
corresponding contour is projected at bottom with a color 
bar which shows intensity of the peaks. Despite a shorter 
distance, the enhancement of intensity at low u for T < Tb is 
weaker than for Pb-Pb G{r,uj) at r = 4 A. 
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